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Abstract 

We propose that general D-dimensional quantum field theories are dual to (D + 1) -dimensional local 
quantum theories which in general include objects with spin two or higher. Using a general prescription, we 
construct a (D + 1) -dimensional theory which is holographically dual to the D-dimensional O(N) vector 
model. From the holographic theory, the phase transition and critical properties of the model in dimensions 
D > 2 are described. 
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I. INTRODUCTION 



Quantum field theory is a universal language that describes long wavelength fluctuations in 
quantum systems made of many degrees of freedom. Although strongly coupled quantum field 
theories commonly arise in nature, it is notoriously difficult to find a systematic way of under- 
standing strongly coupled quantum field theories. 

The anti-de Sitter space/conformal field theory correspondenceyJ43|] opened the door to under- 
stand a class of strongly coupled quantum field theories. According to the duality, certain strongly 
coupled quantum field theories in D dimensions can be mapped into weakly coupled gravita- 
tional theories in (D + 1) dimensions in large N limits. Although the original correspondence 
has been conjectured based on the superstring theory, it is possible that the underlying principle is 
more general and a wider class of quantum field theories can be understood through holographic 
descriptions flfHZO, which may have different UV completion than the string theory. 

In this paper, we provide a prescription to construct holographic theories for general quantum 
field theories. As a demonstration of the method, we explicitly construct a dual theory for the 
D-dimensional O(N) vector model, and reproduce the phase transition and critical behaviors of 
the model using the holographic description. 

The paper is organized in the following way. In Sec. II, we will convey the main idea behind 
the holographic description by constructing a dual theory for a toy model. In Sec. Ill A, using 
the general idea presented in Sec. II, we will explicitly construct a holographic theory dual to the 
D-dimensional O(N) vector model. In Sec. Ill B, we will consider a large iV limit where the 
theory becomes classical for O(N) singlet fields in the bulk. In Sec. Ill C, the phase transition and 
critical properties of the O(N) model will be discussed using the holographic theory. 

II. TOY-MODEL : 0-DIMENSIONAL SCALAR THEORY 

In this section, we will construct a holographic theory for one of the simplest models : 0- 
dimensional scalar theory. In zero dimension, the partition function is given by an ordinary inte- 
gration, 

Z[J] = [ d$e- s[ * ] . (1) 
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We consider an action S[$\ = S M [®] + Sj[$\ with 

S M [$] = M 2 $ 2 , 

oc 
n=l 

Here Sm is the bare action with 'mass' M. Sj is a deformation with sources J n . The values of 
JnS are not necessarily small. In the following, we will consider deformations upto quartic order 
: J n — for n > 4. However, the following discussion can be straightforwardly generalized to 
more general cases. 

For a given set of sources J n , quantum fluctuations are controlled by the bare mass M. One 
useful way of organizing quantum fluctuations is to separate high energy modes and low energy 
modes, and include high energy fluctuations through an effective action for the low energy modes. 
Although there is only one scalar variable in this case, this can be done through the Polchinski's 
renormalization group scheme[8]. First, an auxiliary field $ with mass \i is introduced, 

Z[J) = fi J d^e-^ + ^ 2 \ (3) 

At this stage, $ is a pure auxiliary field without any physical significance. Then, we find a new 
basis <fi and <p 

$ = + 0, 

$ = + (4) 

in such a way that the 'low energy field' has a mass M' which is slightly larger than the original 
mass M. As a result, quantum fluctuations for become slightly smaller than the original field $. 
The missing quantum fluctuations are compensated by the 'high energy field' with mass m . If 
we choose the mass of the low energy field as 

M' 2 = M 2 e 2adz (5) 

with dz being an infinitesimally small parameter and a being a positive constant, we have to 
choose 

MM' m'M 
A = — B = -—, (6) 

m jj, M jj, 



where 



'2 ,,2 e M 

m 2 = M 2 ^— A = -. (7) 

e 2adz _ 1 2adz 



Note that m" 2 is very large, proportional to 1/dz. This is because carries away only infinites- 
imally small quantum fluctuations of the original field $. Moreover, m is independent of the 
arbitrary mass /i because <p is physical. 

In terms of the new variables, the partition function is written as 

Z[J] = + [ AW e-^^+^W^). (8) 



V M' m 



If we rescale the fields, 



<j> ->■ e- adz <P, 4> e~ adz 4>, (9) 

the quadratic action for low energy field <p can be brought into the form which is the same as the 
original bare action, 

Z[J] = mj d<f)d4>e-^ + ^ +M ^ 2+m2 ^\ (10) 

where 

4 

SA<P + 4>] = ^Jn(0 + 0) n , (ID 
n=l 

with j n = J n e~ nadz and m = m e~ adz . Note that j n 's become smaller than the original sources 
J n , which is a manifestation of reduced quantum fluctuations for the low energy field <f>. The new 
action can be expanded in power of the low energy field, 

S, [<f> + $\ = Si [4>] + {ji + 2J20 + 3j 3 2 + 4j 4 3 )0 

+ (j2 + 3j 3 + 6j 4 2 )0 2 + (33 + 4j 4 0)0 3 + J40 4 . (12) 

In the standard renormalization group (RG) procedure [8J, |9J], one integrates out the high energy 
field to obtain an effective action for the low energy field with renormalized coupling constants. 
Here we take an alternative view and interpret the high energy field <\> as fluctuating sources for 
the low energy field. This means that the sources for the low energy field can be regarded as 
dynamical fields instead of fixed coupling constants. To make this more explicit, we decouple the 
high energy field and the low energy field by introducing Hubbard-Stratonovich fields J n and P n , 

Z[J] = m J d<pd^ n=l (dJ n dP n ) e -(s;+AfVW<?) ; (13) 
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where 



+tP 1 J x - iP^jx + 2j 2 + 3j 3 2 + 4j 4 3 ) + M 

+lP 2 J 2 - tP 2 {j2 + 3j 3 + 6j 4 2 ) + M 2 
+lP 3 J 3 - iP 3 (j 3 + 4j 4 0) + J 3 3 

+iP 4 J 4 - iPiji + J 4 4 . (14) 

Now we integrate out <p to obtain an effective action for the source fields. The mass m 2 for the high 
energy field is proportional to 1/dz and only terms that are linear in dz contribute to the effective 
action (for the derivation, see the Appendix A), 

Z[J] = j d^U n=1 (dJ n dP n ) e-mfi+W+smw), (15) 

where 



S w [J, P] = l ( J n ~Jn + nadzJ n )P n 

n=l 

+ + 2P l J 2 + 3P 2 J 3 + AP 3 J 4 ) 2 



(16) 



with J n = (J n +J n )/2. 

After repeating the steps from Eqs. © to (fl~5l ) R times, one obtains a path integral for the 
partition function 

z\j] = [ n« =1 n^ =1 (Dj( fc + 1 ) J DP«) e - sCil) [ J(fe) ' pW ]z[j( ii+1 )], 



(17) 



where 



S<*»[J« P«] = f^i(J^-J^+nadzJ^)P^ 

k=l n=l 

+^(i jf > + 2P« + 3P 2 (fc) jf + 4Pf Jf ) 2 

;(fc) _ / ,(fc+l) , T (k), , r 



(18) 



with Jn = (Jn + J n )/2 and J n = J n - The non-trivial solution for Eq. (flTT) is given by 

= 1 K =1 (DJ n DP n ) e~ s ^ p \ (19) 



where 



dz 



a 



2M 2 



«(5 z J n + naJ n )P n 



[%J X + 2Pi J 2 + 3P 2 J 3 + 4P 3 J 4 ) 



(20) 



Here DJ n DP n represent functional integrations over one dimensional fields J n (z), P n {z) which 
are defined on the semi-infinite line [0, oo). The boundary value of J n (z) is fixed by the coupling 
constants of the original theory, J n (0) = J n . P n (z) is the conjugate field of J n (z). The physical 
meaning of P n becomes clear from the equation of motion for the corresponding source field. By 
taking derivative of Eq. (fl4l) with respect to J n , one obtains 

< (f) n >= -i < p n > . (21) 

Therefore, P n describes physical fluctuations of the operator <p n , and an expectation value of P n 
along the imaginary axis gives an expectation value of the operator. 

The theory given by Eqs. (IT9T ) and (|20l) , which is exactly dual to the original theory, is an one- 
dimensional local quantum theory. The emergent dimension z corresponds to logarithmic energy 
scale[10]. The parameter a determines the rate the energy scale is changed. Despite the apparent 
similarity with the standard RG theory, there is an important difference. In the usual RG approach, 
quantum fluctuations of high energy modes modify coupling constants for low energy modes and 
generated new terms in the effective action. In the present approach, new terms are not generated 
and the structure of the bare theory is maintained at each level. In particular, the highest order 
coupling ( J 4 in this case) obeys the strict constraint, 

<9 2 J 4 + 4aJ 4 = (22) 

which is nothing but the classical scaling. This is because there is no dynamics for the conjugate 
field P 4 for the highest order coupling. Instead, other couplings acquire non-trivial dynamics and 
some of them become propagating modes in the bulk. Fluctuations of those dynamical coupling 
fields embody quantum fluctuations in this approach. 

In quantum field theory, there is a redundancy as to what energy scale one should use to define 
theory. This makes the reparametrization of the RG flow to be a gauge symmetry. Because of 
this redundancy, the partition function in Eq. (fT9l) does not depend on the rate high energy modes 
are eliminated as far as all modes are eventually eliminated. Moreover, at each step of mode 
elimination, one could have chosen a differently. Therefore, a can be regarded as a function 
of z. If one interprets z as 'time', it is natural to identify a(z) as the 'lapse function', that is, 
a(z) = a/ g zz (z), where g zz (z) is the metric. Then one can view Eqs. (fT9)) and ((201) as an one- 
dimensional gravitational theory with matter fields J n . This becomes more clear if we write the 
Lagrangian as 

L = P n d z J n -aH, (23) 
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where H is the Hamiltonian (the reason why H is not Hermitian in this case is that we started from 
the Euclidean field theory). However, there is one important difference from the usual gravitational 
theory. In the Hamiltonian formalism of gravity ll 1111 . the lapse function is a Lagrangian multiplier 
which imposes the constraint H = 0. However, in Eq. (fT9l ), a is not integrated over and the 
Hamiltonian constraint is not imposed. This is due to the presence of the boundary at z = which 
explicitly breaks the reparametrization symmetry. In particular, the 'proper time' from z = to 
z = oo given by 

POO 

1= I a(z)dz (24) 
Jo 

is a quantity of physical significance which measures the total warping factor. To reproduce the 
original partition function in Eq. (OQ) from Eq. (fT9l) , one has to make sure that I = oo to include all 
modes in the infrared limit. Therefore, / should be fixed to be infinite. As a result, one should not 
integrate over all possible a(z) some of which give different I. This is the physical reason why the 
Hamiltonian constraint is not imposed in the present theory. This theory is a gravitational theory 
with the fixed size along the z direction. Nonetheless, the partition function does not depend on a 
specific choice of ct(z) as far as I is fixed. If one wants to make this gauge symmetry more explicit, 
one could integrate in the gauge degree of freedom by summing over different a(z) with fixed I. 
For example, we can integrate over k X1)Z2 parametrizing /-preserving fluctuations of a as 

a(z) = a {z) + k ZuZ2 (5(z - zi) - 5(z - z 2 )), (25) 

where oto(z) is a 'gauge fixed' lapse function. Integration over k Z1:Z2 generates the constraint 

H( Zl ) = H(z 2 ). (26) 

However, this is a trivial constraint which is already implemented in Eq. (fl9l ) with fixed a, because 
it is simply the conservation of 'energy'. 

Although there are many fields in the bulk, i.e. J n , P n for each n, there is only one propagating 
mode, and the remaining fields are non-dynamical in the sense they strictly obey constraints im- 
posed by their conjugate fields. This is not surprising because we started with one dynamical field 
$. There is a freedom in choosing one independent field. In this case, it is convenient to choose 
J3 as an independent field because the conjugate field P3 is multiplied by J 4 in Eq. (|20b . where J 4 
is non-dynamical due to Eq. (1221) . Integrating over P 3 , one obtains 

POO 

S — dz 



M i 
——^(d z J 3 + 3aJ 3 ) 2 - -rrid.Js + 3aJ 3 )(3J 3 P 2 + 2J 2 P X + U x ) 



+ i{d z Ji + aJ x )P x + i{d z J 2 + 2a J 2 )P 2 



(27) 



7 



Pi and P 2 are Lagrangian multipliers which impose the constraints, 

J 2 

(d z Ji + aJi) = —(d z J 3 + 3aJ 3 ), 
(d z J 2 + 2aJ 2 ) = -J.(d z J 3 + 3aJ 3 ). (28) 

4 J4 

Remarkably, these constraints have a local solution, that is, the fields J t and J 2 at a scale z depend 

only on the independent field J 3 at the same scale, 

73 7 

Jl " 16J| 2J? 
3 7 2 

J 2 = V f + J2e~ 2az . (29) 

o J4 

This locality is guaranteed because all source fields at a given scale are tied with one fluctuating 
field <j) at the same scale. Here we considered the case with Z 2 symmetry where J n — for odd 
n. From this one can write down the local action for the independent field J 3 in the bulk, 



/■oo r J./T2 ( T4 rr -2az "| 



(30) 



The first term is a bulk term and the second term is a boundary term. Since J 3 (z = 0) = 0, 
the boundary terms contribute only at the infrared limit z = 00. The theory for J 3 is free in the 
bulk, but the boundary term contains non-trivial interactions. Presumably, this theory is not easier 
to solve than the original theory due to the boundary interactions. However, the construction of 
the dual theory for the toy model illustrates how one can construct dual theories for more general 
field theories. Now, rather than trying to analyze the theory (|30l ), we will move on to apply the 
prescription to more non-trivial field theory : D-dimensional O(N) vector model. 

HI. D-DIMENSIONAL 0(N) VECTOR THEORY 
A. Construction of dual theory 

We consider a D-dimensional vector field theory, 

Z[J\ = J D$ a e -^[*]+^[$]) ; (31) 
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Ja^a + Jab$a$b + Jabc® a$ b® c + Jabcd® a<&b§ c<& d 

(32) 



where 

S M [$\ = y"dxrfy$ a (x)G' M 1 (x-y)$ a (y) 
Sj[$) = f 

+Xld^ a d^ b + Xl^ a d^ b d^ c + xi cd ^ b d^ c d^ d 

Here J dx and J dy are integrations on a D-dimensional manifold Ai D . Here we use Ai D = M 15 
for simplicity. $ a is O(N) vector field. G^-(x) is the regularized kinetic energy with 

G£(x) = J dpp 2 K- l {p/M)e^\ (33) 

where = is an analytic function of s 2 , which remains to be order of 1 for s < 1 

and grows smoothly for s > 1, for example, 

K-!( s ) = e s2 . (34) 

The mass scale M is a UV cut-off above which fluctuations of $ a are suppressed. Sj is a de- 
formation of the free theory. We consider sources J ab ... which are fully symmetric in the flavor 
indices a,b, .... In general, the sources may depend on x. Although we can add more general 
deformations, we will proceed with this quartic action ((32b which is sufficient to illustrate general 
features of the holographic description. 

To integrate out high energy modes, we add an auxiliary vector field $ a , 

Z[J] = [detGs]-^ 2 J D&DQ e-V"W +s *W + ~ s &\ (35) 

where 

S[Q] = y'rfxdy<l a (x)G' D 1 (x-y)<l a (y). (36) 

The form of the propagator Go for the auxiliary field does not affect the final answer. Then, we 
find a new basis <p and 4>, 

$ a (x) = o (x) +0a(x), 

<f> a (x) = j dy (A(x,y)0 a (y) + J B(x,y)0 a (y)), (37) 

where A and B are chosen to satisfy 

G ]s l + A T G D 1 A = G^,, 
G A j + B T G D 1 B — G 1 , 

GJ^+A^B = (38) 



so that the low energy field has a slightly smaller UV cut-off M' = Me adz and the high energy 
field 4> has a propagator G = —(G M > — G M )- Then the partition function can be written as 

Z[J] = [det G'- 1 det G- 1 , det G M ] N/2 [ D(f>D<f> e ~( W^+V M+^D, ( 39 ) 



where 



S' 



A rescaling of x and the fields, 



dxdy a (x)G' 1 (x-y)0 a (y). 



(40) 



x ->■ e aa!2 x, 

a _^ e (2-d)«b/2^ 



->■ e 



(2-d)a<te/2 



(41) 



brings the kinetic energy for the low energy field to the original form as 



where 



with 



Z[J] = [det G~Y /2 J D0D4> 



_ {Sj [4>+4>]+s M [4>}+s[4>]) 



rfx 



ja<Pa + JabMb + jabcMbQc + j abcd4> a4>b4> c4> d 



+3abdi<p a dj<pb + fi c (t>adi<t>bdj(t>c + jl cd <p a <pbdi<p c d. 



(42) 



(43) 



Ja(x) 


= e ^ Qdz Ja(e ad2 x), 




= e 2ad2 J ab (e adz x), 


jabc(x) 


= e^ adz J abc (e adz x), 


Jobcd(x) 


= e^ adz J abcd (e adz x 






•7aL( X ) 


= ^ adz Jl(e adz *), 


iafecd( X ) 


= e^ adz J% cd (e adz x) 



(44) 



and 



G-\x - y) = e (2+d)ad2 G''- 1 (e ad2 (x - y)). 
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(45) 



The new action can be expanded in power of the low energy field, 



Si 



+ 



+ 

+ 
+ 
+ 



dx{ |j a + 2j ab (j) b + 3j abc (f) b (f) c + 4:j abcd (f) b (f) c (f) d 

- 2a i (%d j $ b ) +&Wc - 24(jIWc) 

+ Zjabcd&di&cdjfa - 2d i (j l J bcd 4> b 4) c d j 4> d ) (p a 

jab + 3j abc c + Qjabcd<t>c<t>d ~ d^f^dj^c) 
+ fabcd^cdj^d ~ Idiijabcd^cdrfd) 



>a9b 



Jab + jabc^c + Jabcd<Pc(pd 



di<t) a dj4>b 



jabc + ^jabcd<Pd ~ ^ditiabcddj^d) 



>aVbVc 



jabc + ^jabcd'Pd 



^adi^bdj^c 



(46) 



+jabcd<Pa<Pb<Pc<Pd + jabcd^bdi^cdj^dY 

Here we ignored boundary terms in the D-dimensional space, assuming that the boundary is at 
infinity where couplings vanish. We decouple the low energy field and the high energy field by 
introducing Hubbard-Stratonovich fields, 

Z[J] = [detG- 1 ]^ 2 / D</>D$DJDP e-rt +s "M +§ l*», 



(47) 



where 



= sm 

J d,X {iP a J a ~ iP a 



ja + 2j ab (j) b + 3j abc 4> b 4> c + 4:j abcd (f) b (f) c (j) d 



+ Ja<Pa 



+iPab-Jab — iPab 



jab + 3j abc (f) c + Qjabcd<Pc<Pd ~ di(j a i c dj4> c ) 



+ faicd d i4>A^d - 2d i (j a J bcd 4> c d j 4> d ) 



+ Jab<t>a4>b 



+iPab,ijJ a i ~ iP ab,ij j ab + Kbc^c + j a \ cd <t>c<t>d 



^riPabcJ abc iPi 



abc 



jabc + ^jabcd<Pd - ^di(j l a J bcd dj(f) d ) 



+ J- b Ui(p a Uj(p b 

+ Jabc4>a4>b 



+iPabc,ijJabc ~ iP abc,ij [jabc + 2 jabcd<t>d\ + f^a^b^c 
^P abed J abed — iPabcdjabcd + J abcd4> a<Pb<\> c<P d 
+iPabcd,ijJ a bcd ~ iPabcd,ij3abcd + Jabcd^bdi&cdjfa. 



(48) 
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Now we integrate out <fi, following the similar procedure as explained in the Appendix A to 
obtain 



Z[J] 



d(j)DJDP e 



-(5jM+SmM+S (1) [J,P]) 



(49) 



where 

S m [J,P] =Jdx 



lPa{dJ a 



+iP abc (dJ t 



2 + D 



aj a ) + iP a b(dJab ~ ^OiJ ab ) + iPab,ij(dJ%) 



abc 



Q — D {j 2 — d 

-aJabc) + ^Pabc^[p J abc 7^^J C 



abc' 



+ \jd*dy |s a (x)G(x-y)s a (y)|. 
Here, d = £ - a J2f =1 Xi £- and 



-iPabcd{dJ a bcd - (4 - D)aJ abcd ) + iPabcd,ij{df a \ cd - (2 - d)aj^ cd )\dz 

(50) 



iJa + ^PbJab — ZdjiJ'abdiPb) + ^PbcJabc ~ ^ji^JnPiPl 



be 



+ Pbc,ijJaL + ^PbcdJabcd - -^djiJ^diPbcd) + ^J % abcdPbcd,ij 



(51) 



with J a _, = (J a .. + J a .)/2 and Tj = (J£ + J l a {)/2. If one integrate out J's and P's in Eq. 
( |49l , one reproduces the action obtained by integrating out <p in Eq. (|42l) to the order of dz. If we 
keep applying the same procedure to the action for the low energy field as we did in the previous 
section, the partition function can be written as 



Z[J] 



DJDP e - s[J > p] 



(52) 



where the (D + 1) -dimensional action is given by 

S[J, P] = J dxdz \iP a {dJ a - ^-^-aJ a ) + iP ab {dJ ab - 2aJ ab ) + iP^idJ^) 



6 - D 

+iP a bc(d.J abc - aJabc) + iPabc,ij(d.J a { c 



2 - d 

-^- aJ abc) 



+iPabcd{dJ a bcd - (4 - D)aJ abcd ) + iP a bcd,ij{dJ a { cd - (2 - d)a.J a 3 bcd ) 



+ \jdyidydz |as a (x)G'(x - y)s a (y)j. 



(53) 



with 



iJ a + 2P b J ab — 2dj(J a 3 b diP b ) + 3Pb c Jabc ~ 9j(J a J bc diPbc 

2 
3 



+ Pbc,ijJ a bc + ^PbcdJabcd ~ q ^ (^L^i-Pfccd) + 2J a lcdPbcd,ij 



(54) 
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and G (x) = M3mGm{^)- Here the partition function is given by the functional integrals of the 
source fields J and their conjugate fields P in the (D + 1) -dimensional space Ai D x [0, oo). If 
the P-dimensional manifold Ai D has a finite volume V, the volume at scale z is given by Ve~ aDz 
because of the rescaling of space in Eq. (|4T|) . Since Gj\/(p) _1 is smooth in momentum space, 
the last term in Eq. ([53]) can be expressed in gradient expansion. Higher derivative terms are 
suppressed by (f|) n . Therefore the full theory is local in (D + 1) -dimensional space. The locality 
along the emergent z direction is due to the fact that physics at a given energy scale E depends on 
higher energy physics only through physics at an infmitesimaily higher energy scale Ee adz . 

As is the case in the 0-dimensional theory discussed in the previous section, a can be regarded 
as the lapse function in a gravitational theory in (D+l) dimensions, but the Hamiltonian constraint 
is not imposed because of the boundary at z = 0. Moreover, there is no fluctuating shift function 
which imposes the momentum constraint Pj = in the usual gravitation theory. This is, again, due 
to fact that the boundary explicitly breaks the diffeomorphism symmetry and only the momentum 
conservation (not p = 0) is implemented in the theory. 

One key difference from the 0-dimensional theory is that there exist bulk fields with non-trivial 
spins. In Eq. (1531 , there are spin two fields which are coupled to the energy momentum tensor at 
the boundary. In the presence of more general deformations in the boundary theory, one needs to 



introduce fields with higher spins II 12l 



B. Large N limit 

In this section, we will consider a large N limit where the dual theory becomes classical for 
O(N) singlet fields. To see this, one decomposes tensors with rank two or greater than two into 
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traceless tensors with the same rank and tensors with lower ranks as 

Jab = J2<)ab + -jyJab, 

1 1 - 

Jabc = [Jzahc + J^ac + J^c^ab] + J^Jabc-, 

Jabcd = -j^(5ab5 c d + 5 ac 5bd + 5 a d5bc), 

Pab = -P^afc + Pab, 

Pabc = P^bc + P^ac + P^ab + -Pa&c, 

Pabcd = P^ab^cd + ^ac^bd + ^ad^bc), 

jJb = J^dab+^jJb' 

I 1 - 

^aL = [^3a^6c + ^35 #ac + J 3c ^ab\ + JpJ a bd 

i ■ J 4 
Jabcd = -J-j-($ab8 c d + 8 ac 8bd + dad^bc), 

Pab,ij P2,ij&ab Pab,ij 1 

Pabc,ij Psa,ij$bc P^b,ij^ac Pzc,ij^ab Pabc,iji 

Pabcd,ij P^,ij 

{Sabred + 8 ac 5 bd + 5 ad 5 bc ), (55) 

where the fields with the bar are traceless, J aa = J aab = J l J a = J l J ab = and P aa = P aab = 
P aa ,ij = P aa b,ij = 0. For the quartic couplings which are the highest order couplings, only the 
O(N) invariant parts are kept, assuming that the bare quartic couplings are O(N) invariant. Note 
that the structures of the highest order coupling fields are not modified at all z, which is not true 
for other coupling fields : for lower order couplings, non-singlet contributions are generated at low 
energy scales even though the bare couplings are O(N) invariant. Now, we consider the large N 
limit with fixed J a , J 2 , J ab , J 3a , J abc , J 4 , J 2 J , J*, Jg, J a { c , J 4 , f 4 3 . In this limit, the leading order 
action becomes 

S[J, P] = J dxdz [iP a {dJ a - ^-^-aJ a ) + iNP 2 {dJ 2 - 2aJ 2 ) + i^P ab (dJ ab - 2aJ ab ) 

+iNP 2iij (dJ i 2 j ) +i±P abjij (dJ a j b ) 

6 D 1 6 D 

+i3P 3a (dJ 3a — aJ 3a ) + ij-^Pabc{dJabc ^ — ocJ abc ) 

+i3P 3a ^{dJl - ^ajg) + i^Pabc^idJll - ^aJ a l) 
+i3NP A (dJ A - (4 - D)aJ A ) + i3NP 4jij (dJ l 4 j - (2 - d)aJ l 4 j )} 
+ lf dxdydz {as (x, z)G'(x - y)s (y, 2)}, (56) 
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where 



lJ a + 2P a J 2 + ^P b J ab - 2dj{jjdiP a ) - ^djifidiP,) 

6 3 

+3P 2 J 3a + —PabJsb + -J^PbcJabc 

-d 3 (jgd t p 2 ) - ^-(JS&Aa) - ^- 2 dA J lAPb C ) 

■ 2 1 

I p f!J I _p ffl j p p3 

r 2,ij J 3a jy- r ab,ij<J3b ' jy 2 bc M J abc 

+12P 3a J 4 - 28,(^8^) + (ir.'!':^ 



(57) 



In the large N limit, the action is manifestly proportional to N. Therefore, one can ignore 
quantum fluctuations for singlet fields such as J 2 and J 2 , and it is enough to consider saddle 
point solutions for singlet fields to compute correlation functions of singlet operators. It would be 
natural to integrate out all non-singlet fields and obtain an effective theory for single fields alone. 
However, it turns out that the effective action for single fields become non-local in this O(N) 
vector model. This is because there are light non-singlet fields in the bulk and integrating over 
those soft modes generates non-local correlations for singlet fields. This means that we should 
keep light non-singlet fields as 'low energy degrees of freedom' in the bulk description if we want 
to use a local description. 



C. Phase Transition and Critical Behaviors 



In this section, we will describe the phase transition and the critical properties of the model 
in D > 2 using the holographic theory. To maintain the locality of the description, we will keep 
light non-singlet fields in the theory by choosing J3 a 's as independent fields. We will focus on 
the simple case where there is no deformation on the energy-momentum tensor, and all sources 
respect the O(N) symmetry, 

J a (x,0) = 0, 

J afe (x,0) = J 2 (x)$ab, 

J abc (x,0) = 0, 

J abcrf (x, 0) = ^ {SabScd + S ac Sbd + SadSbc) (58) 

with J l 2 j = j l i h = Jg = J l J bc = fl = 0. From now on, it will be assumed that J 4 (x) = J 4 > 
is ^-independent, but J7~ 2 (x), which may have either sign, is rr-dependent in general. We first 
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integrate over P 3a to obtain 



S[J, P] = I dxdz [iP a {dJ a - a Jo) + iNP 2 {dJ 2 - 2aJ 2 ) + i-^P ab {dJ ab - 2aJ ab ) 

s 6 — D 1- - 6 — D - 
- ijj-{dJ 3a 2 — aJ ^ + i j^ p abc(dJ abc — aJ abc ) 



+ i3NP 4 {dJ 4 - (4 - D)aJ 4 )j 
+^ J dxdydz C(x,2;)[aG'(x-y)]- 1 C(y,z)}, 

where C(x, 2) = (<9J 3a (x, z) - ^aJ 3a (x, s)) and 



2 6 3 - 

iJ a + 2P a J 2 + —P b J ab + 3P2^3a + -TrPahJsh + -TT^PhcJahc 



(59) 



(60) 



Now we integrate over P a , P 2 , P ab , P abc , P A to obtain constraints, 



dJ a - 2+ 2 D aJ^j - -^j- (^dJ 3b - 2 ~° aJ 3b ^j (^2J 2 5 ab + -^J ab 

N(dJ 2 - 2aJ 2 ) - — (dJ 3a - -^—aJ 3a j J 3a 
3 A_ 6-£> \ 

C^3a " a^3a ^36 



4iVJ 4 



N \ ah ah) 
OJib 7, 0£j 3 b J 3a 



4N 2 J 4 



0J 3 a - 



2 

6- D 



0. 
0. 



OtJ 3c Jahc = 0, 



6 - D - 

dJahc CtJahc = 0, 

<9J 4 - (4 - D)aJ A = 0. (61) 



For the boundary condition (1581) . the solution for the constraints is given by 



J a (x, z 



J 2 {x,z 

Jab (X, Z 
Jabcip^-i Z 

J 4 (x, Z 



3 w , w , , , J 3 a(x^)J 2 (e-x)e 2 - 

J 3a (x, Zj J 3 b(x, Zj J 3f >(x, 2) H - f ~ , 



16iVJ 4 (x,z) 2 

1 ^3a(x, z) J 3a (x, z) + (e QZ x)e 2az 



2 J 4 (x, z) 



87VJ 4 (x,z 
3 



4 J 4 (x, z) 
0, 



J 3a (x, z)J 36 (x, z) 



(62) 
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FIG. 1: Saddle point configuration of Js a {z) (a) in the disordered phase and (b) in the ordered phase. 
Although the action is quadratic in J 3a in the bulk, the boundary action drives the phase transition. When 
J 2 is sufficiently negative, a Mexican-hat potential at the boundary drags Jz a (z) away from J$ a {z) = in 
the bulk. At the critical point, J 3a at the IR boundary z = oo is more or less free to fluctuate, generating 
algebraic correlations between fields inserted at the UV boundary z = 0. 

where all repeated indices are summed as usual. Using this result, one can eliminate all fields 
except for J 3a to obtain the action, 



In the first line of the above expression, all fields which do not have an explicit argument are 
understood to be at (x, z). The first term is a boundary action and the second term is the bulk 
action. In the bulk, J 3a is massless and its action is purely quadratic. Although the bulk theory 
is non-interacting, the boundary term has non-trivial interactions and quantum fluctuations for the 
non-singlet field are still important. 

Due to the boundary term, the theory has a non-trivial phase transition. In the boundary term, 
the only contribution is from the infrared boundary, zi = oo because there is no symmetry break- 
ing source in the bare theory, J 3a (x, 0) = 0. To discuss the phase transition, it is enough to 

consider x independent singlet source fields J 2 and J A . If J 2 > J 2 C , where J 2 C J A M^ D ~^ 

is a critical value, the boundary potential has the minimum at J 3a = as is shown in Fig. ED (a). 




(63) 
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Therefore, the minimum energy configuration of J 3a is the straight line along the ^-direction. On 
the other hand, if Ji < J% the boundary term has the Mexican-hat potential. This causes the bulk 
field J 3a to deviate from the trivial configuration J ia = 0, leading to the spontaneous symmetry 
breaking jl^j as is shown in Fig. [H(b). This describes the second order phase transition from the 
disordered phase to the ordered phase. We can use —iP^a ~ 2 0a as an order parameter which 
roughly measures the slope of J 3a (;z) along the z direction. In the phase with broken symmetry, 
J 3a in the bulk is proportional to yN and fluctuations away from the saddle point is 1/N sup- 
pressed, leading to the mean-field like critical exponent for the order parameter, (3 = 1/2 in the 
large iV limit. 

Now let us compute the two-point correlation function of the 0(N) singlet operator <p 2 at 
the critical point. For this we integrate over J 3a in the bulk in the presence of ^-dependent source 
J2 (x) = J~2+J^ 2 ( :x ) with the boundary condition J 3a (x,0) = 0. Since the bulk action is quadratic 
one can use the saddle point solution in the bulk. In terms of the Fourier modes, 



Jpxe az + (6-D)az/2 



(64) 



J 3a (x, z) = J dp f a (p,z)e 

the bulk action becomes 

Sbulk = J dzd P e- 2az [G']- l (pe az )[d z fa(p,z)][d z f a (-p,z)]. (65) 

The equation of motion for / a (p, z) is 

d z (e- 2az G'-\pe az )d z Up, z)) = 0. (66) 

There are two independent solutions : 1 and K(pe az /M). The boundary condition at z = fixes 
only one coefficient, leading to one-parameter family of solutions given by 

K{pe az /M)~ 



f a (p,z) = y a (p) 



(67) 



K (p/M) 

Here K(s) is the regulator defined in Eq. (|34l) . and y a (p) is an arbitrary function. This arbitrariness 
is due to the fact that we have not imposed the boundary condition at the IR limit z = 00. The 
condition of finiteness of the action alone does not fix the IR boundary condition because both 
independent solutions are regular in the IR limit. Here one should integrate over y a because y a 
is also a dynamical field. In this way, the boundary condition in the IR limit is dynamically 
determined by the boundary action, whereas the boundary condition in the UV is fixed by the bare 
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coupling. Using the saddle point solution, one obtains an action for the boundary field y a , 
1 



i6Ji 



16J 4 



dp p 2 K'\p/M)\y a (p)\ 2 

^2 J dp!dp 2 Jftfa + p 2 ) + ^(-Pl - Pa) ya{Vl)ya{V2) 



+ 256iY\7 3 / rf P lrf P2 rf P3 rf P4 ^(Pl + P2 + P3 + P4)2/a(pi)2/a(p2)ya(P3)Z/a(p4)- (68) 

Remarkably, this action has the same form as the original action for <f> a , although this action is 
for the source field at the IR boundary. Despite the fact that it is no more convenient to use this 
dual description than using the conventional field theoretic method for the O(N) model, one can 
proceed to integrate out the boundary field y a , by summing over the usual chain diagrams, to 
obtain an effective potential W^fq)] f° r me singlet field. From this, one can readily compute 
the correlation function. For example, in D = 3 one obtains 

d 2 W 

W7P~, V^Tp-f V ~ q + (analytic terms) (69) 

This leads to the scaling dimension [<p 2 ] = 2 as expected. 

The fact that one has to solve essentially the same field theory to dynamically impose the IR 
boundary condition in the dual description is somewhat disappointing from the perspective of dual- 
ity. However, this may have been expected in the 0(N) model. This is because all N independent 
fields are massless, and we kept all of them in the local bulk effective action. Since no dynamical 
field has been integrated out, no information of the original theory has been coarse-grained and 
the bulk theory carries the exactly same amount of dynamical information as the original theory. 
An alternative approach would be to include more singlet fields such as J 2n associated with the 
operators (4> 2 ) n for n = 1, 2, N, and integrate out all non-singlet fields to obtain an effective 
action for iV independent singlet fields. However, if one keeps only a small number of singlet 
fields and integrate out the remaining fields, the resulting action will not be local because their 



masses in the bulk are more or less equally spaced ! 15H . 

The holographic description for the O(N) model is not very useful for practical purpose. After 
all, we can solve it in a much simpler way. We believe that the holographic approach may become 
more useful for more strongly interacting theories where there exist a gap between a small set of 
operators with small scaling dimensions and a large set of operators with large scaling dimensions. 
In such models, it is expected that quantum fluctuations of heavy propagating modes in the bulk 
will be encoded in local effective actions for light modes which carry dynamical information only 
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for operators with small scaling dimensions. Then, imposing the IR boundary condition for the 
light modes may become more tractable than solving the original field theory. 



IV. SUMMARY AND DISCUSSION 



In the present paper, we provided a general prescription to construct holographic theory for 
general quantum field theory. Through an explicit construction, we showed that the holographic 
description for the O(N) vector model correctly reproduce the spontaneous symmetry breaking 
phase transition and critical behaviors, as predicted by standard field theory methods. In the future, 
it is of great interest to apply this method to other strongly cou pled quantum field theories where 
standard field theory techniques fail, such as matrix models ll 1 6M 1 811 and non-Fermi liquids in 2+1 
dimensions which are expected to be dual to certain matrix models 111 911 . 

If general quantum field theories in D dimensions are holographically dual to certain (D + 1)- 
dimensional local theories, the next question would be "What (D + 1) -dimensional theories are 
dual to local D-dimensional quantum field theories ?" If one knows the answer to this question, 
one may use the holographic description to define quantum field theories and to classify them. To 
answer this question, it may be helpful to fully understand general structures behind the present 
holographic construction. 
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VI. APPENDIX A 

The action in Eq. (fl~4~l) can be expanded in power of the high energy field as 

4 

S'j = Sj[(j)] + i( J n - j n )P n 
n=l 

+0'i - 2zPij 2 - 3iP 2 j 3 - UPdS 
+(j 2 - ZiPija - 6*P 2 j 4 )0 2 
•(./:; • l//';./.;)^ 

+J40 4 



Integrating over to the order of 1/m 2 , one obtains 



v. = / 



m 



d<J)Tl 4 n=1 (dJndPn) _ 

V ra l + B 



- (m V+gj M+iELi V»-jn)P»+ 4( J 2 +B) ) 
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where 

A = (ij x + 2P 1 j 2 + 3P 2 j 3 + 4P 3 J 4 ), 

B = (j 2 -3iPij 3 -6iP 2 u). (72) 

The cubic and higher order terms in do not contribute to the linear order in 1/m 2 ~ dz. If we 
keep only those terms that are linear in dz in Eq. (1711) , we obtain the action, 

8 = M 2 2 + SM + z £(J n - j n )P n + ^- 2 + ^~r (73) 

n=l 

However, it is not easy to take the continuum limit (dz — >■ 0) in this expression for the following 
reason. We can regard j n and J n as being defined at coordinates z and z + dz respectively, where 
z is the logarithmic energy scale in the renormalization group flow. Then P n is defined in the 
interval (or at the middle point of the interval), [z, z + dz}. Usually, A 2 dz can be interpreted as 
the integration of A 2 between z and z + dz in the continuum limit. This would be the correct if 
A were a fixed constant of the order of 1. In the present case, however, A contains the dynamical 
field P„ whose typical amplitude is order of m ~ h=. Therefore, an error of order of 1/m in j n in 
the integrand gives a non-trivial contribution to the integration, leading to a discrepancy between 
the result in Eq. (1731) and the one obtained in the naive continuum limit. 

To fix this problem, we consider the following trick. First we absorb the factor ( - m2 1 fB ) m front 
of A 2 in the action into the measure of P3; we change the variable P 3 to P3 in Eq. (fTTT) . 



which leads to 



where 



m 4mj4 



Z[J] = y#n n ^(rfJ^Pn)^3rfP3e-( M ^ 2+SjM+5 '), (75) 



A' 2 

n^3 



. / \/m 2 + B , y/m 2 + B — m , . 

+i(Js ~ js) 1 ^3 + , • ^ + 2P ^ 2 + 3P2 ^ 3 (76) 

\ m 4mj 4 1 

with A' = (iji + 2Pij 2 + 3P2J3 + 4P3J4). If we take the leading order terms, the above expression 
becomes 

A '2 

S' = I J2(Jn ~ 3n)Pn + ^ + ^ ~ h)^ (77) 
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Dropping the prime sign in P 3 , the partition function becomes 

Z[J) = J d<pii n=l (dJ n dP n ) e~ s \ (78) 

where 

4 i 

S" = M 2 2 + 3j[<f>] + l ( J n - 3n)Pn + ^{iji + 2Pij 2 + 3P 2 j 3 + 4P 3 j 4 ) 2 . (79) 

n=l 

However, this expression is not completely satisfactory either. If one integrates out J n and P n in 
this expression, one obtains an action for <p which is different from the result one obtains after 



integrating out (p directly from Eq. (fT2l) to the order of em boli . The difference is the contribution 
from the tadpole diagram. The tadpole diagram simply shifts the local couplings, and it turns out 
that its contribution can be accounted for by replacing j n with (j n + J n )/2 in the last term of Eq. 
(EU) as 

4 i 

S'" = M 2 4> 2 + SM + *( J n - Jn)Pn + ^(ui + 2Pj 2 + 3P 2 j 3 + 4P 3 J 4 ) 2 , (80) 

n=l 

where j n = [j n + J n )/2. Although this is an infinitesimal change, it still gives a non-trivial 
contribution because (J n — j n ) ~ 0(l/m) and P n {J n — j n ) ~ 0(1). It is straightforward to 
show that with this action one reproduces the same action as the one obtained by integrating out <fi 
directly from Eq. (fT2l) to the order of dz. In Eq. (l80l) . the mean value of j n and J n is multiplied to 
P„ . Just as the error of the trapezoidal method in the usual discrete integration is suppressed to dz 3 , 
the difference between Eq. (f8~0l) and the integration in the continuum limit becomes sub-leading in 
dz even though P n ~ 1/ \fdz. Therefore, Eq. (|8~0l) can be readily extended to the continuum limit. 
If we use j n = J n e~ nadz and keep the leading order term, we obtain Eqs. (fT5l) and (fT6l) . It is noted 
that if one uses Eq. d73l and take the naive continuum limit, some couplings are spuriously shifted 
due to the amplified error introduced in the continuum limit. 
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